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The ■paper expands the basic concepts of the coding theorists in the 
representation oj data strings by algebraic polynomials, and develops the 
representation oj both time and location of individual binary positions in 
the same polynomial. Further, it advances a set of algebraic operations on 
such polynomials to correspond to the various subfunctions that are ac- 
complished in the actual domain circuits. The specific applications of the 
techniques proposed in this paper for the design and synthesis of such circuits 
is presented in a companion paper. 

"Choose a set of symbols, endow them with certain properties and postulate 
certain relationships between them. Next, . . . deduce further relationships between 
them .... We can apply this theory if we know the "exact physical significance" of 
the symbols. ... The applied mathematician always has the problem of deciding 
what is the exact physical significance of the symbols. // this is known, then at any 
stage in the theory we know the physical significance of our theorems. But the 
weakest link of physical significance is extremely fragile." The original source of this 
principium is J. E. Kerrick in An Experimental Introduction to the Theory of Pro- 
bability, Belgisk Import Company, Copenhagen. It is also quoted in a slightly 
different form by F. M. Heza in An Introduction to Information Theory, McGraw- 
Hill Book Co., New York, 1061. 

I. INTRODUCTION 

Magnetic domains exist freely in thin platelets of orthoferrite crystals 
obtained by slicing them so that their crystalline axis is perpendicular 
to the surface of the platelets. Such domains are also present in very thin 
epitaxial garnet films (Fig. la) on suitable substrates. When the platelets 
or films are subjected to bias fields, these domains assume cylindrical 
shape and their diameter shrinks to microscopic sizes (Fig. lb). Such 
domains (also called "bubbles") are stable under an appropriate bias 
field condition and they may be manipulated to perform' ' 2 storage, 
gating, looping and also certain elementary logic functions. 3 

The domains are generally propagated from one location in the circuit 
to the next by subjecting them to the local bias field gradient. Basically 

1535 



1536 



THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1972 



there are two methods of providing such a field gradient to propagate 
the bubbles. In the "field access propagation," 4 an alternating magnet- 
ization is imposed in a patterned soft magnetic overlay by an in-plane 
rotating magnetic field generated by a pair of coils carrying an alter- 
nating current. The coils completely surround the platelet with their 
axis in its plane. Two of the most commonly used overlay patterns are 
shown in Figs. 2a and b. During one cycle of the alternating current in 
the coils, the domains in the platelet move from one point in a pattern 
to the corresponding point in the adjoining pattern. This finite distance 
that the domain traverses during one cycle is defined as a "period." In 




(a) 




(W 



Fig. 1— (a) Magnetic domains as they are observed by Faraday effect in a typical 
epitaxial film 5 to 8 microns deep, deposited on Gadolinium-Gallium-Garnet (GGG) 
substrate 20 to 40 mils thick. Magnification 340. (b) Formation of "bubbles" from 
magnetic domains at a bias field of 30 ()e in same material used in Fig. la. Magnifica- 
tion 340. 
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Fig. 2 — (a) T-bar type of overlay used for field access propagation, (b) Y-bar type 
of overlay used for field access propagation. 



the "conductor propagation, " the local field gradient to move the 
bubbles is supplied by a current in a conductor (Fig. 3). A single phase 
current is generally used to periodically shift a bubble from one stable 
position to the next. Such stable positions are derived by a soft magnetic 
overlay also embedded on the platelet, and the periodicity of movement 
of the bubble from one position to the next depends on the frequency of 
the single phase excitation of the conductor. The distance which the 
bubble moves during one cycle of the single phase current is also defined 
as one "period." 

Typical orthoferrites (YbFeO ; , , YFe0 3 , etc.) sustain 40 to 50 micron 
diameter bubbles, and the period is approximately 200 microns. Typical 
garnets (E^T^Ah^Fe.^Oiz and Gd2.:,Tbo.7Fe30 12 ) can support 4 to 8 
micron diameter bubbles and the period is about 25 microns. The 
orthoferrites require about one micro-second to shift a bubble position 
by one period. The newer garnet materials also require about the same 
time, thus yielding a data rate of about one megacycle. It is customary to 
employ "bubble-no-bubble coding" with field access propagation and 
"lateral displacement coding" (LDC) with conductor propagation. In 
the former type of coding, the presence or absence of a bubble at an 
appropriate location denotes one or zero. In the latter type of coding, 
the bubble positions are coded as one or zero by laterally displacing 
them from one coding position to the other coding position (see Fig. 3). 

All the bits of information are propagated by one period in one clock 
cycle in the field access propagation. In conductor drive circuits with 
lateral displacement coding, the average velocity during propagation is 
generally limited to one finite value, even though information bits are 
sometimes held stationary. When one finite velocity of propagation is 
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Fig. 3 — Conductor propagation of bubbles. 



assumed, it is possible to extend the capabilities of the conventional 
polynomial algebra used by coding theorists 50 to encompass the space- 
time relationships in bubble circuits. Bubble circuits perform in the time 
dimension and in the space dimension. The space dimension is further 
divided into subdimensions; the different sections (or elements) of the 
circuit which perform independently. Hence, any algebraic representa- 
tion should encompass the representation of time and the representation 
of space dimensions which constitute the circuit. It is proposed that the 
time dimension be associated with X and the space dimension be 
associated with Y in the algebra. 

II. DECOMPOSITION OF TIME, SPACE (CIRCUITS) AND FUNCTIONS 

Consider the subclassification of time and space dimensions as follows: 

The total time for a circuit to perform a function consists of a series of 
individual time intervals neccessary for the submodular functions. 
These individual time intervals can each be represented as a certain 
known number of clock cycles. Thus the unit of time is one clock cycle 
at the excitation frequency of the main field in field access propagation, 
or is one clock cycle at the drive circuit frequency in conductor 
propagation. 

The physical layout of the circuit can be classified into various sections 
(elements such as paths, loops, functional modules, etc.). Each element 
further consists of a certain predefined number of periods. This leads to 
the unit of physical (or spatial) dimension as one period corresponding 
to one pole pitch in the T-bar, Y-bar or chevron pattern in field access 
propagation, or to one pole pitch of the driving conductor in conductor 
propagation. 
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Next consider the space-time relation. In field access propagation, 
the speed of the bubbles is one period per clock cycle, and is influenced by 
the angular velocity of the main rotating field; and stationary bubbles 
are rarely encountered.* In conductor propagation, the frequency of the 
exciting circuit determines the velocity;* and stationary bubbles are 
commonly encountered. However, in every circuit, over a limited 
duration and within a preselected element of the circuit, one can express 
the space-time relation with absolute certainty. 

Finally, consider the overall algebraic representation. The entire 
circuit function is modeled by a series of algebraic operations, each one of 
which corresponds to a subfunction in the circuit. Each subfunction is 
carried out in the time dimension and in the space dimension. Hence, if 
we can resolve the function into its subfunctions, the circuit into its 
elements, and time into sets of clock cycles, and identify the individual 
subfunction with the circuit element and the appropriate set of clock 
cycles, then we can analyze and predict the functioning of the circuit 
with great accuracy. The representations of individual subfunctions by 
corresponding algebraic operations are developed in Section IV. 

III. REPRESENTATION OF TIME, LOCATIONS AND BINARY VALUES OF A BIT 
POSITION 

The origin of time may be chosen to be at any desired instant. How- 
ever, a certain amount of flexibility and ease of representation results if 
the origin of time is chosen to coincide with a definite function in the 
circuit. Generally, bubble circuits perform repetitive functions and it 
may be convenient to choose the origin of time at the start of a repetitive 
cycle. When circuits perform a wide variety of nonrepetitive functions, 
then the analysis should be attempted for each function independently 
to ascertain the correct operation of each one of the functions. In the 
algebraic analysis of bubble circuits, it is proposed that the exponent 
of A' (the carrier of time dimension as introduced earlier) be used to 
represent the number of clock cycles that have elapsed between a 
prechosen origin of time and the instant under consideration. 

Further, it is proposed that the location of any given binary bubble 
position be represented by two components: (/) the element in which the 
binary position is presently located and (ii) the exact period in that 



* The momentary variations of bubble speed at crossovers and compressors are 
ignored, and the entire distance is considered as one period. The effect of corners is 
dealt with separately. 
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element at which the binary position is currently present. This leads to 
designating a subscript to Y (the carrier of space dimension as intro- 
duced earlier) to indicate the element number and an exponent of Y to 
indicate the exact location within that element. 

Finally, it is proposed that the binary value of a bubble position at a 
given instant of time be denoted by a. The binary value changes as a 
bubble position passes through the known transition points in the circuit. 
However, at a given time, which is a certain number of clock cycles past a 
preselected origin of time (a known exponent of X), and at a given 
location (known subscript and known exponent of Y), the binary value of 
a bubble position is cither known, or it can be determined with absolute 
certainty from other circuit considerations. 

3.1 Representation of an Isolated Bubble 

Examine a single binary bubble position (Fig. 4a) the binary value of 
which is 'a' at an instant of time j clock cycles from a prechosen origin of 
time within the fcth element of the circuit located at the Zth location. 
Then it may be represented as 

u = aX'Yl. 



1 2 3 



12 3 
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Fig. 4 — (a) Represent a (ion of a single bubble. See text for explanation of j, k, 
and /. (b) Representation of a bubble stream j clock cycles after a prechosen origin of 
time. 
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3.2 Representation of a Bubble Stream 

Consider a string of n binary bubble positions (Fig. 4b). Let the binary 
values of these positions be a„, a, , ■ ■ • , a„_, during a clock cycle which is j 
clock cycles from a prechoscn origin or time. If these binary positions 
are located in the fcth element at location /„,/,,-••, /„_, , respectively, 
then the string of data can be represented as 

u = OnX'Fi' + • • • a^.X'Fl."- = X i "JT ai Y['. (1) 

1=0 

The sign of individual terms does not carry any significance. When the 
binary positions are adjacent to one another, then /„,/,, • • ■ , /„_, are 
consecutive numbers. 

Example 1: Consider four bubble positions the binary value of which 
is a ,a^,a 2 and a 3 . If they occupy the 3, 2, 1, and location of a sixth 
element after 28 cycles from a prechosen origin of time, then they may be 
represented as 

u = X 2 *(a ( ,Yl + • • • a a Yl) = X 28 £ a,r<-°. (2) 

i=0 

3.3 Explanation oj the Algebraic Representation 

From the point of view of circuit analysis, the algebraic representation 
leads to the following propositions. 

(i) A series of bubble streams in a circuit are represented by a series of 
polynomials. 

(ii) Any one bubble stream is represented by a particular polynomial 
(the sum of individual terms). 

(Hi) Each bit within a bubble stream is represented by a term (the 
product of components). 

(iv) The binary value of the bubble position is represented by a. 

(v) The number of clock cycles between a prechosen origin of time and 
the end of the cycle under consideration is the exponent of A". 

(vi) The location of the bit of information is represented in two 
sections: the element within the circuit (the subscript of F) and the 
location within the element (the exponent of )'). 

3.4 Implications of Representation 

The implications of the prechosen representation are: 

(i) That a "snapshot" (i.e., a complete description of binary values of 
bubble positions and their respective locations) may be extracted from 
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the algebra after a predetermined interval of time from a prechosen 
origin of time (a known exponent of X). 

(ii) That there is a unique polynomial for every bubble stream in a 
circuit during any one cycle. 

(Hi) That all the pertinent information regarding all binary bits of 
information in a stream is available within the algebraic polynomial 
representing it. 

IV. REPRESENTATION' OF FUNCTIONS 

4.1 Generation of Bubble Streams 

Choose an origin of time at the end of the generation cycle of the first 
bubble position (i.e., as it is leaving the generator). If the generator is 
going to generate n binary positions the values of which are a , a, , ■ • • , 
and a B _, , then at the end of (n — 1) clock cycles, the binary string just 
generated may be written as 

u = X-^aoY*' 1 + ■ • • a n ^Y°) = X"~ l "j? a.-F"" 1 " (3a) 

t'-O 

at an instant when the a,,-, position is just leaving the generator. The 
alternate representation of the string after n cycles is 

u = X-'Z a,F ( "-°. (3b) 

i-O 

In (3a) and (3b), the exponents of Y indicate the locations along the 
bubble path which lead out of the generator. When the coefficients 
flo , a, , • • • , a„_i are consistently one, the action of an unconditional 
generator is represented (see Fig. 4 of Ref. 4, and Fig. 3 of Ref. 3 for 
T-bar and Y-bar configurations). 

4.2 Annihilation of Bubble Streams in Bubble-No-Bubble Coding 

The function of annihilation of any bubble stream u may be simply 
represented as 

u = 0. 

It is important to note that it is not identical to a polynomial in which 
all the binary bit positions are zero. When a conditional generator 
generates a string of binary bits which are all zero, it is still necessary to 
represent the binary string (3a) or (3b), since this string of data may 
interact with other strings at a later point in the circuit. 

The function of a conditional annihilator may be represented as the 
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action of a conditional gate (to be discussed later) and the action of an 
unconditional annihilator. 

4.2.1 Resetting of Streams in Lateral Displacement Coding 

In a functional sense, the function of resetting in Lateral Displacement 
Coding (LDC) is analogous to the function of annihilating in bubble-no- 
bubble coding. However, there is an important representational dif- 
ference. The location of a bubble position which is reset in LDC is still 
identifiable, whereas the bubble position in field access circuits loses its 
identity and location upon entering the annihilator. Upon resetting in 
LDC, the bubble positions enter a "dead interval", or a sort of coma 
from the activity of the circuit, yet the algebra has to account for the 
elapsed time during which the bubble positions are in the reset state. 
Hence, a stream of n reset bubble positions in locations 0, 1, 2, • • • , 
1 may be represented in the usual way as 



n 



u - X* Z a,}"" 1 -'', (4) 

t'-O 

where ; represents the number of clock cycles from a prechosen origin of 
time, and where a n , a t , ■ • ■ , a„_, invariably represent the reset status. 
Example 2: The representation of 3 reset bubble positions in location 9, 
5 and 1 after 20 clock cycles past a given time origin is 

u - X> F 9 + aX° + a 2 Y l ) = X 20 £ a<F (9 - 4 ° . (4a) 

4.3 Functions of One Stream Resulting in One Stream 

Let u p be the initial polynomial about to undergo the function, and let 
u q be the resulting polynomial after m clock cycles. In a categorical sense, 
the operation may be represented as 

U Q = F m (u p ); or, F m (u p ) -> u„ , 

where F may represent: temporary freezing, translation in the forward 
direction, translation in the reverse direction, looping, special looping 
used in memory operations for dynamic data allocation, etc. 

4.3.1 Temporary Freezing of Bubble Streams in their Locations 

This function is generally encountered in conductor pattern and rail 
propagation with lateral displacement coding. Consider an n bit data 
stream ;' clock cycles from a prechosen origin of time, and represented as 

u p = X i '' jf aY l k , (5) 
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where a and I are a, and U respectively. After freezing the movement of 
the stream for m clock cycles, we have the stream represented as 

u a = X m -u p = X i+m 'j: aY' k . (6) 

If the bubble positions are all located at adjoining locations then 

1=0 

where (I) is the location l of the bubble position a . 

4.3.2 Movement o] Bubble Streams 

This is the most common bubble function and it is necessary to estab- 
lish a sense of directionality in the movement. If the bubble stream is 
moving so that a n - x bubble position of eq. (1) occupies position l after a 
certain number of clock cycles, then a„ would be the leading bubble. It is 
easier to work with a„ as the leading bubble,* and this implies that 
lo > h > h ■ • ■ > l*-i ■ If the positions are in adjoining locations, then 
Z,-_, is l t + 1. Consider a stream of bubbles represented by u v in eq. (5), 
which has moved in the forward direction for m clock cycles, then the 
resulting polynomial u q is 



t-n-1 



u a = X m Y m -u p = X i+m £ aY' k +m , (8) 

i=0 

where a and I as defined earlier are a, and I, ; . If the data positions are in 
adjoining locations with a in (I), then 

u Q = X i+m 'T! aYl l) ' i+n . (9) 

t' = 

Generally, the bubble stream crosses elemental boundaries when it 
moves. If this stream shifts from element k to element t during the 
movement, and if their intersection is located at Y' k and Y° t (Fig. 5), 
then 

u a = X i+m '"E *YT-* (10) 



* This notation helps the circuit designer to comprehend the location and move- 
ment of the leading bubble position first, rather than comprehending the location and 
movement of the last bubble position first. 
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Fig. 5— Forward movement with a shift from element k to element /. 

in general.* If the bubble positions are located in adjoining positions 
I — h = Iq — i = (I) — i and 



u Q = X i + m Z aY[ l) - i+m -. 



(11) 



Example 8: The representation of 7 bubble positions, the binary 
positions of which are a through a fi located at the 15th, 14th, ■ • ■ , 9th 
locations of the 3rd element of a circuit at an instant of time, 36 clock 
cycles from a prechosen origin of time, is 



= x 36 Ea.-rr 1 '. 



(10a) 



Further, if this stream travels from element 3 to element 4, after travers- 
ing for 20 cycles with the boundary between elements 3 and 4 being 
located at Y\° and F.° , then the final polynomial is represented as 



but 



t=0 
F 35-.- = j.35-,-20 = yu- t 



* An alternate way to visualize the crossing of boundaries is to write 
< = A'' + "> "t 1 aY k ' +m , 

i-0 

and then replace n»+- by Y ,", yielding a = / + m - z, which leads to eq. (10). The 
prime indicates that the polynomial as such does not represent an observable stream 
but will do so after the next operation (s). 
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and thus 

i = 

The exponent of F 4 becomes negative if the stream does not traverse at 
least 11 cycles when it is 7 bits long, or if the stream traverses for 20 
cycles but is 16 bits long. It is thus possible to have the stream segmented 
as u v and u Q by an improper choice of the traversing time. Sometimes 
negative exponents may be carried for a few steps in the analysis without 
interpreting polynomials as streams during these steps. 

Examine a bubble stream represented by u p in eq. (5) having moved 
in the backward direction for m clock cycles, then 

u Q = X m Y- m -u p = X i+m 'xTa..^-, 

t=0 

where a and I are o, , and h respectively, and a a is the last bubble in this 
case. If the bubble stream moved from element k and entirely shifted 
in tth. element, then 

u q = X i+mi Z^Yr l - m) , (12) 

1=0 

where* the intersection of elements k and t is located at Y° k and Y\ . If 
the bubble positions are in the adjoining locations, then l t = Z — i 
where Z is the location of the bubble a in u p before its transformation 
to u Q . 

Example 4: The representation of a 4 bubble stream, the binary values of 
which are a , a v , a 2 and a 3 located in the 7, 5, 3, 1 periods of element 5 
at an instant 35 clock cycles after a prechosen origin of time, is 

u v = x™ E«.^r 2i - ( 12a > 

i = 

After 9 cycles of backward movement, if the stream is in element 3 
with the boundary of elements 5 and 3 located at 7° and Yf , then 

u Q = X" £a f 75 8 -". ( 12b ) 



* An alternate way to visualize the crossing of the boundary is to write 

< = X^ m "if 1 aiY ,'-">, 

and then replace Y k ° = Y,' +a , when a is < 0, thus yielding z + a = z + (I - m) 
which leads to eq. (12). 
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4.3.3 Looping of Bubble Streams 

Looping is generally encountered in T-bar or Y-bar propagation when 
a fixed delay or storage is necessary for bubble streams without actually 
freezing their movement. 

Consider a bubble stream represented by u„ in (5), which is circulating 
in a loop with a periods* for m clock cycles, then 

u a = X m Y m -u„ = X i+m ' °£ aFl + ", 

i = 

but in a loop (/ + m = (I + "0 mod g) thus leading to 

u Q = X i+m "2 an i+m)aode , (13) 

i-0 

where a and I are read as a, and /, respectively, and /, = (/„ — i) if the 

bubble positions are in the adjoining location with a as the leading 

bubble. 

Example 5: The representation of 26 bubble positions 1 a„ , a, , • • • , a 25 , 

occupying locations 28, 27, • • • , 3 in a loop (element 5) at an instant 

1214 clock cycles from a prechosen origin of time is 

u P = X 1214 Ea.-rf--. (13a) 

i=U 

After looping for 980 clock cycles, the final representation is 

u Q = X 2194 '£ a.-yr 08-0 mo ' 130 = X 2194 2 a.T? 3 -''. (13b) 

i=0 i=0 

4.3.4 Movement Around Corners 

Corners in T-bar and Y-bar circuits need special attention since 
bubble positions lose or gain a quarter period when a 90-degree turn is 
present. In most cases, their effect may be eliminated by considering a 
movement to span two or an even number of compensating 90-degree 
corners. However, when it is necessary to predict the movement of 

* The number of T-bar periods should be considered as the number of dock cycles 
to bring back the leading bubble to its original location in the loop. The orientation 
of the T-bars and the direction of rotation both play an important part in the deter- 
mination of the number. In any case, the value of g is the actual number of T-bars ±1 
depending on the orientation of T-bars and the direction of rotation. 

f One encounters this bubble stream in (39, 26) the shortened BCH encoder with 
magnetic domains constructed along the same principles as discussed in Section 3 of 
Ref. 7. 

* When the effect, of the corners in a loop is being considered, the value of k used 
should be (h mod 4), and the number of periods in the loop should be measured as 
indicated in Sec. 4.3.3. Else the value of A- should be its real value. 
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streams which include k number of 90-degree unidirectional corners, 
then a polynomial u„ in eq. (5) becomes 

Uq = rr" M -m, (14) 

where the sign is determined by the direction of the turn with respect to 
the direction of rotation of the main driving field and the orientation of 

the T-bars. 

No special algebraic consideration is necessary in conductor pattern 
and rail propagation. 
Example 6: Consider the two bubble positions represented as 

u p = X 20 (a n Yl + o,H). (14a) 

If these positions traverse for 13 cycles, and encounter three unidirec- 
tional -90-degree turns in the path, then the final polynomial is 

u a = X 33 y (,3 -* ) -u p = X 33 (a o rr l + a, Yi 8J ). (14b) 

The fractional exponent of Y indicates that the binary positions a and 
a, are lagging 270-degrees behind their corresponding positions had 
there been no turns. If it is necessary to bring them back in phase, then 
an additional £ clock cycle is required, and the polynomial u Q would be 

u a = X 33i (a Yl° + a,Y 1 :). (14c) 

4.3.5 Routing oj Bubble Streams 

This function plays a critical role when it is necessary to transfer 
streams of binary information into a certain branch of a circuit at a node. 
The algebraic representation after this function is identical to the 
polynomial u a in eq. (10), when the polynomial crosses the boundary of 
one element k at Y' k and enters another element t at Y° t . In the general 
case, when the bubble stream is channeled from an element p at Y' p , 
and enters another element t at Fi" , w„ may be represented as 

u q = X^'lf a,y! + — +< ' > , do) 

i=0 

where a and I as defined earlier are a, and /, respectively. 

4.3.6 Inverting the Binary Content of Bubble Streams 
Consider a polynomial 

u P = Z*'"E a tp Y l h . (16) 

i=0 

If the bubble stream so represented has gone through an inverting gate 
during the following m cycles, then the resulting polynomial may be 






BUBBLE CIRCUITS 1549 

written as 

m, = X 1 *" "e" a iq Y' k +m , (16a) 

1=0 

whore / denotes U and 

a„ = d, p . 

Example 7: Consider 3 data positions in an LDC circuit. If a n , a t , a 2 
denote their status at an instant 49 clock cycles from a time origin, and 
are located in the seventh element at locations 17, 13 and 9 then, 

u p = X 40 I>„,rr 4 '. (16b) 

i = 

If these binary positions pass through an inverting gate located at F™ 
and Y\ , then the bubble stream after 15 clock cycles is 

u„ = X ni 2>*«r""" (16c) 

1=0 

4.3.7 Opening and Closing Gaps in Bubble Streams 

This function is quite effectively used in dynamic data reallocation 
with T-bars. The bubble stream in the loop has two preferred paths, one 
for each direction of rotation. In one direction, the stream traverses 
(n + 1) periods in the loop, and in the other, it traverses (n) periods in 
the loop. (See Fig. 6). If after z clock cycles of clockwise movement 
of a bubble stream in the loop the direction is reversed for one clock 
cycle, then the data bit at the nth (Fig. 6) location is at the Oth location, 
and the data position which was at (0 — 1) mod n) location is at the «th 
location. Now z clock cycles of anticlockwise rotation would have 
effectively included the data position at the ??th bit at the (z mod (n + 1)) 
location, and all the remaining positions one location behind their 
original locations. A converse process takes place for excluding a data 
bit position in the stream. 

The algebraic equivalent of this function may be represented as 
follows: 

Let the numbering of the locations in the loop be in the direction 
(anticlockwise in Fig. G) which permits the stream to traverse n periods 
in the loop. Let the contents of the loop (without the nth position) be 
represented as 

Mp = A" "if ai Y\ (17) 
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Fig. 6 — Schematic representation of the dynamic data reallocation function. 

The polynomial u* undergoes a shift in clockwise direction for m (with 
m < n — 1) clock cycles. The resulting polynomial u q becomes 



i = 

As a distinct step from the above representation, let 



(18) 



(19) 



represent the contents of the loop including the nth location before a 
series of counter clockwise shifts for m clock cycles. The resulting poly- 
nomial is 



= X i+m J^ ai Y (i+ ' a) 



(20) 



Example 8: Consider a string of bubbles the binary values of which 
are a , a, , ■ ■ ■ , a 7 at an instant 92 clock cycles from a time origin 
located at 0, 1, 2, • • • , 7 in a memory loop (Fig. 6) with 255 (i.e., n — 1) 

* There is no need for a subscript for Y since only one element (the loop) is being 
considered. 
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periods in the clockwise direction and 256 periods in the anticlockwise 
direction. This can be represented as 

Up = X 92 f>,r. (17a) 

1=0 

Let the 256th (i.e., nth) position contain the bubble position a n at the 
same instant of time. Four cycles of rotation of the field causing a 
clockwise circulation of binary positions yields 

I i II 1=4 ) 

In the first four terms in the polynomial (38), (i — 4) mod 256 would 
correspond to 252, 253, 254 and 255 for i ranging from to 3. 

Now if the field is rotated for 5 clock cycles in the opposite direction, 
resulting in an anticlockwise shift of bubble positions, then the resulting 
polynomial is 

u = Z ' 0, {i: a,F (2S7+,,mo ' 1257 + a n Y 4 + £ a,T' + 1 } , 

I i=0 i=4 ) 

or 

u = X 101 {£S a.F + a„7 4 + £ a.y , + 1 )- (20a) 

Vi=0 i=4 J 

In effect it is seen that the bubble position after the nth location has been 
inserted at the 4th location in the bubble stream. 

4.4 Functions of One Stream Resulting in Two or More Streams 

Let u p be the initial polynomial about to undergo the function resulting 
in two polynomials u„ and u, . In general, functions of this type may be 
represented as 

u Q + u T = F m (u p ); or, F m (u p ) -» u Q + u r , 

where F* represents duplication or addressing sections of the initial 
polynomial u„ into one or the other branch of a circuit. A finite number 
of clock cycles (m) are allowed during which the operation takes place. 

4.4.1 Duplication and Replication of Bubble Streams 

Duplicators may require a certain finite number of clock cycles to 



* These functions are used extensively in the companion paper 8 dealing with the 
applications of the algebra. 
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Fig. 7 — The bubble stream u p entering the duplicator D. 

operate. In the T-bar propagation, one clock cycle is necessary to dupli- 
cate. The effect may easily be included in the algebra by considering the 
duplicator as an extra period of element p with u p (see eq. 5), and the 
zero location of elements q and r in which u Q and u T will be positioned 
after m clock cycles. (See Fig. 7). 

Consider a string of binary data represented as u p in eq. (5) to be 
duplicated during the m clock cycles yielding u Q and u r represented as 



u„ = X Y -Up = X 2-j aY « 

u r = X m Y m - d -u p = X i+m "if aY l r +m — d , 



(21) 
(22) 



where a and / are a, and h respectively, the duplicator is located after z 
periods in k, and d is the number of clock cycles to accomplish the 
duplication. It is to be noted that when the bubbles occupy adjoining 
locations, d cannot exceed 1 for satisfactory duplication. Under such 
conditions* 



-i+m-'-l 



u i — x 2-j aY " 

i = 

and a similar expression for u r , where (I) denotes l . 



(23) 



* An alternate way to visualize this transformation is to consider that the dupli- 
cator is located at the intersection of Y p ' +l , Y q ° and Y r ° which leads to Y p " = 
Y q <*-u+u = Y r a ~ i,+U , when a > z + 1, thus yielding 



u' = X'' 



"£ l aVyn-.+m, 



where (I) is the location l of the first bubble position a prior to duplication. This 
leads to u„ in eq. (23). 
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Replication leads to a series of resulting polynomials u„ , u r , ■ ■ ■ , etc., 
but the algebraic treatment is exactly the same. 

4.4.2 Addressing Sections oj a Stream into Different Branches of a Circuit 

Consider a polynomial u p in element p which will approach a gate. 
The gate operating for g clock cycles will address the first n' bubble 
positions in q and the rest into r. This function can be visualized as the 
effect of two independent translatory functions: (see Section 4.3.1) (i) 
the first n' positions are translated in the forward direction and change 
elements from p to q, and (ii) the remaining positions are translated in 
the forward direction and change element from p to r. Algebraically this 
can be expressed as follows: 

= X**£ aY l p = V V aV ! : |- A" /\ n\\, (2-0 

with 

(25) 



and 



"£ aY l p = X j ' if aY' p + X* 

i = 1=0 


i-n' 


t=n'-l 

vit" V" 1 „V ,+m -' 

u Q = X 2^i a *Q 

i=0 




i=n-l 

u r = X i + m 2 aY***-', 





(26) 



where a and Z are a, and l< respectively (see Section 5.3.1). The gate 
should be located at Y' p , Y° and Y° T . It can be seen that the first bubble 
position does not reach the gate till (z — (/<> + 1)) clock cycle, and to 
divert the first n' positions the gate should be operating to divert into 
element q for exactly (l Q — l n - + 1) clock cycle, leading to the design 
detail that the gate should divert into q for g clock cycles where g = 
(z — l n - + 1). If the bubble positions are in the adjoining location, then 
the gate has to operate diverting into q for n' clock cycles starting after 
(z — (l Q + 1)) clock cycles. Further, it has to act for the next (n — n') 
clock cycles to divert the bubble position into r. 

When a gate addresses various sections of a data stream into more 
than two elements, the algebraic representation is similar. Such a con- 
dition exists if m is not chosen large enough in the previous case, and 
then there will be two bubble strings u„ and u r together with a section 
of u p , which has not been processed by the gate. 
Example 9: Consider a data stream* in element of a circuit. Eight data 

* Such a data stream is encountered in general rate change circuits represented 
in Fig. 1 and la of Ref. 9. 
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bits a through a 7 , 12 clock cycles after a prechosen origin of time, 
occupy location Y x2 through Yl . The polynomial representing the stream 

is 

u p = X 12 g aH 1 "'. (24a) 

1 = 

If the data stream passes through a gate at Y i3 which diverts the first 4 
bits into element F", and the last four into Y 2 , then after 14 clock cycles 
the resulting polynomial u q in 1 and u r in 2 may be represented as 



and 



X 26 £a,F! 3 -', 



l T = X 26 2 a «+4^2 '• 



(25a) 
(26a) 



4.5 Functions oj Two or More Streams Resulting in One Stream 

Consider two streams u p and u„ (Fig. 8) interacting to yield one 
stream u T at a gate g), where (|) may denote the function of logical 
gating, combining, etc. This function is different from the previous 
functions, since the individual binary values a, in the polynomials are 
likely to be changed by this function. Let 

u v = X* 'if a iv Y l p (i *\ (27) 



and 



= X< £ a iq Y': i9 \ 



(28) 
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Fig. 8 — Gating or merging of u p and u q to yield u r . 
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where l(ip) and l(iq) represent /,„ , and l iQ respectively. The number of 
locations including the first and the last binary position of u„ is (/„„ — 
£ (n -i. P ) + 1); (i.e., n if adjoining locations). It can be seen (Fig. 8) that 
the length of the resulting bubble stream, u r , is the largest number n" 
of the following four numbers* 

(i) Up ~ lu-i.v) + 1, 

(»") ka ~ iu'-i.,) + 1, 

(m) Qob ~ !<»-..*>) + 1 + 2 - K - (z' - /,„) 

= (z- /,„_,,,) - (z' - U + 1, 
or 

(iv) (/„, - /<„■-,.„) + 1 + z' - Zo. - (2 - U 

= (Z' - ?<,'-,„)) - (2 - /n„) + 1, 

where 2 and z' denote the number or periods along elements p, and q at 
which the gate (g) is located (see Fig. 8). The polynomial u r can be 
written as 

u T = X i+m<mn jt l a ir Y l T lir) , (29) 

where l(ir) represents l, r respectively, and 

o.> - <*.„ (D a,, (30) 

Ur = li V + m - z, (31a) 

or 

l ir = l iQ + m - z'. (31b) 

The subscripts for a and / should be chosen with adequate care to 
consider only the interacting binary positions in streams u„ and u Q that 
pass through the gate (g) simultaneously. It is important to note that 
each term in u r results from a term in u v , or in u Q , or from terms in both. 
When there is no term in one polynomial (u q or u p ) corresponding to a 
particular term in the other (u„ or u Q ), then the appropriate equation 
(31a or b) for the exponent of }' should be chosen. When there is a term 
in one polynomial (m, or u„), corresponding to a given term in (u p or u Q ), 

* (i) or (/'/) indicates the length u,, or u q with the longer stream u p or u q completely 
overlapping u„ or u„ respectively. (Hi) or (iv) indicate partial overlap, with u q or u p 
being nearer the gate. 
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then the location calculated from l iv or /,-, will yield the same result for 

h, ■ 

Example 10: Consider a bubble stream u„ with four binary bits a 0p , 
a lp , a 2p and a 3p , after 41 clock cycles from an origin of time located at 
3, 2, 1 and of element 5 and represented as 

u p = X 4 ' i>^rr\ (27a) 

i = 

Consider a next stream u q with seven bubble positions u„ q through u Bq , 
also 41 clock cycles from the prechosen origin of time occupying 10, 9, 
S, • - - , 4th locations of element 6, and represented as 

u Q = X il £««.*••"'. (28a) 

i=0 

If u„ and u q approach an exclusive-or gate (g) located at Y* 5 , Y l a l , and 
Y° 7 , then the polynomial u T after 12 cycles can be derived as follows: 
Oo, , Oj, pass through the gate during the first three cycles while u p is 
still approaching the gate. The binary positions a 2q through a 5q interact 
with a 0p through a 3p in the gate (g) for the next four cycles. a 6q passes 
through the gate and the gating is now complete. During the last five 
cycles, the bubble stream in element Y 7 moves away from the gate thus 
leading to the resulting polynomial 

u r = X M Za ir Y 7 l -\ (29a) 

i = 

where a ()r , a, r and a nr are a 0q , a iq and a liq respectively, and 

a ir = a, p © a iq (i - 2 through 5). (30a) 

This example corresponds to the case (?) in Section 4.5. 

4.6 Functions oj Two or More Streams Resulting in Two or More Streams 

This function, though rarely encountered in normal bubble circuits, 
can still be conveniently represented as an integral procedure of many 
subfunctions in which two or more streams result in one stream. If g } , 
(72 , f/a • • • are individual functions yielding streams 1, 2, 3 • • • , etc., 
then the algebraic representation of Section 4.5 may be extended to 
represent streams 1, 2, 3 • • • , etc. One such example is presented in 
Section 2 of Ref. 8. 

V. OVERALL BUBBLE CIRCUIT FUNCTIONS 

We have a set of mathematical tools to predict the binary values and 
locations of individual bubble positions as the binary streams undergo 
different submodular functions within the circuit. The interval of time 
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chosen for these submodular functions in Section IV is 'm' clock cycles. 
To study the overall circuit function, the function is divided into a series 
of submodular functions (F, , F 2 , F a ■ ■ ■); the circuit which performs the 
functions is divided into a series of elements 

(1,2,3 ••• ,*,*+ 1, ••■), 

and the time necessary to accomplish the function is divided into a series 
of (m, ,m 2 ,m 3 ■ • ■ , etc.) clock cycles. 

5.1 Subdivision o] Circuit into Elements 

After isolating the subfunctions within the overall circuit function, 
the elements that accomplish these subfunctions may be identified. A 
series of fine functional subdivisions may be necessary to identify the 
particular circuit elements. 

Some of the examples of the elements are storage paths, transmission 
paths, loops, etc. Sometimes a particular element (k) of a circuit designed 
is very short, and it cannot accommodate the entire string of data. When 
it is still desired to study the contents of the n bit binary string of the 
polynomial u p in that element, then it is possible to fictitiously extend the 
element to just accommodate the n data bits. With the observer located 
at a preselected period (say Y b k ) in the element, the binary values of data 
which flow past this location would still be the values of a„ , a x , a 2 ■ ■ ■ a n 
in the calculated polynomial u,. . The instant of incidence of the leading 
bubble a,, would be (n — b) clock cycles prior to its value as predicted by 
the exponent of X(i.e., j n ) in the polynomial u k . 

5.2 Subdivision of Time 

Subfunctions are accomplished by elements within specified intervals 
of time. The interval of time for a specific subfunction is almost entirely 
determined by the circuit parameters and the clock frequency. Any one 
particular interval of time may, however, be conveniently expressed as a 
certain number of clock cycles. 

These different values of clock cycles m, , m 2 , m 3 , • • • , etc., are 
necessary to calculate the polynomials u x , u 2 , u A , • • • , etc., and they in 
turn uniquely define the values, intervals and positions of binary bits in 
the circuit. Generally, the subfunctions may proceed in series or in 
parallel, and different bubble streams may simultaneously undergo 
different functions in different elements of the circuit. However, the 
complete function of a circuit starts and finishes at an instant of time. 
Hence, whenever subfunctions are accomplished in parallel, it should be 
realized that the summation of w, does not equal the total number of 
clock cycles necessary for the complete circuit function. Each sub- 
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function proceeding simultaneously should be modelled independently 
by its corresponding algebraic operation. 

VI. CONCLUSIONS 

The operation of bubble circuits may be effectively analyzed by multi- 
dimensional polynomial algebra without actually constructing the 
circuits. The location of all the data positions in the circuit can be 
accurately predicted at any preselected instant of time during the 
operation of a circuit by this technique. When all the circuit parameters 
are not accurately known, the analysis helps in the calculation of some 
of the circuit parameters. Further, it helps to algebraically check the 
validity and effectiveness of a conceived circuit in the performance of 
specified functions. 

The algebra may also be used for circuits that do not perform instantly 
but need a certain predetermined duration for movement, duplication, 
gating, etc. Several technologies (magnetic domain, charge coupled and 
charge transfer technologies) presently being developed in the Bell 
System fall into this category. 
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